Effective lagrangian approach to structure of selected nuclei far from stability 
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The convergence of the effective field theory (EFT) approach of Furnstahl, Serot and Tang to 
the nuclear many-body problem is studied by applying it to selected doubly-magic nuclei far from 
stability. An independently developed code, which can incorporate various levels of approximation 
of the chiral effective lagrangian, is used to solve the self-consistent relativistic Hartree equations. 
Results are obtained for ground-state properties such as binding energies, single-particle level struc- 
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ture, and densities of the selected spherical, doubly-magic nuclei "■'"'"""Sn and ""''"Ni. Calculated 
spectra of neighboring nuclei diff'ering by one particle or one hole show agreement with the most 
recent experimental data. Predictions for nucleon densities are presented. 
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I. INTRODUCTION 

The study of nuclei far from the /3-stability line has 
drawn the attention of nuclear physicists. The implica- 
tions of the structure and reaction cross-sections from 
these nuclei play an important role in the description 
of nucleosynthesis processes such as the r-process, which 
takes place in the region of neutron rich nuclei. In addi- 
tion to this, these nuclei are the perfect place to test the 
validity of some successful methods and theories, which 
have rendered satisfactory results along the stability line 
for medium-heavy nuclei. The development of new ex- 
perimental facilities and the implementation of new tech- 
niques have already given experimental data for some of 
these, like the doubly-magic nuclei ^5QSn82 and ^ggSnso. 
Some of the new approaches being tested are those of 
relativistic mean-field models. In this paper we are go- 
ing to deal with a relativistic mean-field theory derived 
from the chiral effective lagrangian of Furnstahl, Serot 
and Tang §. This la grangian, which incorporates the 
symmetries of QCD, particularly a non-linear realization 
of chiral symmetry, is built as an effective field theory 
and contains the lowest lying hadronic degrees of free- 
dom as the main ingredients. Though it reproduces the 
ground-state properties of stable doubly-magic nuclei, its 
predictability far from the /3-stability region has not been 
tested, especially the convergence of the calculations as 
the number of effective terms in the lagrangian increases. 

Relativistic mean-field models can account for a large 
number of bulk properties of nuclei and there are dif- 
ferent methods that have produced acceptable results. 
Recently, Toki, et al. ^ used a relativistic mean-field 
approach to calculate properties of nuclei from the pro- 
ton to the neutron drip lines, extending the calculations 
to super-heavy nuclei. The good agreement with known 
experimental data shows how relativistic mean-field mod- 
els can account for a variety of phenomena in the nuclear 
many-body system. Toki and his co-workers, using a 
mean effective lagrangian, have an extensive program un- 
derway on applications of this approach to astrophysics. 
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in particular to element formation in supernovas |3| . 

The idea of using effective field theories to describe 
the nuclear problem is an increasing field of study, an ac- 
count of recent results can be found in Q . The relativistic 
mean-field approach to the nuclear many-body problem 
has been successful in describing medium-size nuclei, and 
its success can be understood under the perspective of ef- 
fective field theory (EFT) and density functional theory 
(DFT) H, ^. The success of previous relativistic mean- 
field models, like the original Walecka model H, and its 
extensions to include nonlinear couplings of the a field, 
is now understood by applying the ideas of naturalness 
from EFT Q . A comparison of various mean- field models 
and an account of their success was presented by Furn- 
stahl, et al 1^. In that paper it was shown, using the EFT 
approach and NDA, why various mean- field approaches 
have been successful and that the use of different degrees 
of freedom to describe nuclear properties corresponds to 
different organizational principles. Thus a description 
of the nuclear many-body system based on nucleon and 
meson fields, as in this approach, or an expansion in nu- 
cleon densities, as in the Skyrme model, are equivalent. 
The principles of EFT rely on the fact that there exist 
natural scales and in order to understand a particular 
phenomenon it is enough to use probes that can resolve 
the dynamics at that scale. These probes are the relevant 
degrees of freedom needed to describe the system. The 
most general lagrangian density is built from these inter- 
acting fields paying attention to symmetries imposed by 
an underlying theory (QCD in this case) and arranged 
by some organizing principle, like the naive dimensional 
analysis (NDA) from Georgi and Manohar [|o[ 0. The 
unknown dynamics, of heavier degrees of freedom, is in- 
tegrated out and appears as coupling constants in the 
theory. These are determined by fitting the theory to 
known experimental data. The number of terms in this 
EFT lagrangian is infinite and, in order to give predictive 
power to the model, it is necessary to have an organizing 
principle that groups together terms that contribute at 
the same level. In this way it is possible to truncate the 
expansion at some arbitrary order and be certain that 
all non-redundant terms, of the same order, have already 
been taken into account. 

The inclusion of new interaction terms in the la- 
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grangian aims to give a better description of the en- 
ergy functional, thus going beyond the Hartree approx- 
imation and incorporating exchange-correlation effects. 
The new terms are non-renormalizable, a distinction with 
standard RMF models, which truncate the effective la- 
grangian at this level without any underlying physical 
justification. A chiral effective lagrangian has been de- 
rived by Furnstahl, Serot, and Tang ||]|, which incorpo- 
rates a non-linear realization of chirial symmetry and uses 
the lowest lying hadronic levels: i) pions, whose field ap- 
pears as the phase of a chiral rotation of the identity ma- 
trix, ii) nucleons, described by isospinor field N, iii) an 
isoscalar- vector field, whose mass is taken to be equal to 
the w-meson mass and iv) an isoscalar-scalar field which 
simulates correlated two-pion exchange with an effective 
mass around 500 MeV, and v) an isovector field, the p- 
meson, to take into account the symmetry energy in nu- 
clei. By applying the relativistic mean-field approxima- 
tion, a set of equations of motion can be derived from 
the effective lagrangian in which the meson fields act as 
classical fields. 

The connection to DFT is made by considering the 
expansion of the relativistic mean-field lagrangian as an 
expansion of the energy functional of the nuclear many- 
body problem in terms of nuclcon densities and auxil- 
iary classical meson fields. According to the principles of 
DFT ||fj, ^ , knowledge of the exact energy functional 
can be used, by optimizing it with respect to the fermion 
densities and auxiliary potentials, to calculate the ex- 
act ground-state scalar and vector densities, energy, and 
chemical potential for the fully interacting many-fermion 
system These optimal observables are used here 

to compare the experimental values with the calculated 
ones. Other observables are not expected to be well re- 
produced, but a comparison with some of them is in- 
cluded in this paper in order to understand the extent to 
which the theory can be applied, as well as to find new 
directions in which more work needs to be done. In this 
way the EFT principles combined with DFT methods 
give an approach to the nuclear many-body problem that 
uses the simplicity of solving the self-consistent Hartree 
equations in which new contributions coming from the in- 
teracting fields can be incorporated in a systematic and 
controlled way. The expansion parameters that make 
the truncation of the effective lagrangian possible in the 
nuclear structure case are essentially the ratios of mean 
isoscalar, scalar and vector fields to the nuclcon mass 
(w 1/3) and the ratio of the Fermi momentum to the 
nucleon mass (« 1/4). 

In this paper we study the convergence of the chiral ef- 
fective lagrangian (GEL) approach to the nuclear many- 
body problem when extrapolated to selected doubly- 
magic nuclei far from the /^-stability region. The study 
focuses at how well additional levels of approximation of 
the energy functional reproduce the above ground-state 
properties of the selected, spherical, doubly-magic nuclei. 
In addition to this, we also examine the effect of fitting a 
larger set of observables, in the mean, along the stability 



line. We center our attention on the four doubly-magic 
nuclei ^|QSn82, ^ggSnso, 28^^20 and 2iNi5o chosen for their 
relevance to the nucleosynthesis processes. Here we com- 
pare the predictions of the GEL model at various levels of 
approximation with the most recent experimental results 
for these nuclei. We compare results of nuclear binding 
energies for these nuclei and also extend this analysis to 
neighboring nuclei that differ by one particle or one hole 
from the doubly-magic nuclei thus comparing the values 
of the chemical potential. We find in both cases remark- 
able agreement with recent experimental data. In addi- 
tion to this we look at the convergence for single-particle 
levels, both for neutrons and protons. We examine how 
sensitive these calculated spectra are to the various levels 
of approximation in the GEL, as new interaction terms 
come into play. Nuclear densities are also shown for the 
"proton rich" nucleus ^5oSn5o and the neutron rich one 
^|QSn82. These densities are in very good agreement with 
other calculations using the relativistic Hartree Bo- 
goliubov equations. According to the EFT approach used 
here, the systematic inclusion of interaction terms will 
approach the true energy functional, going beyond the 
Hartree level and including exchange-correlation terms. 
Thus, here we include additional pairing interactions only 
to the extent they are systematically included through 
the additional terms in the GEL and the fitting of the 
constants. The goal of the present work is to study how 
well the current level of approximation in the GEL repro- 
duces the ground-state properties of the selected nuclei 
without explicit use of the pairing interactions. 

The calculations are done using an independently de- 
veloped code, which solves the self-consistent Hartree 
equations at any level of approximation. The code has 
been tested successfully by reproducing previous calcula- 
tions. The coupling constants, for each level of approxi- 
mation, were obtained by Furnstahl, et al. ||l| by fitting 
in the mean to properties of stable doubly-magic nuclei 
leading to different parameter sets according to the level 
of approximation ||l|, ^ . We use here these results and use 
the same name for the parameter sets, i.e. Wl, Ql, Q2, 
Gl and G2 in increasing level of approximation, respec- 
tively. In addition to these, we include two additional 
sets: L2 and NLG. The parameter set L2 corresponds to 
the same level of approximation as the set Wl, it con- 
tains no non- linear terms in the meson fields. The other 
parameter set, NLG, has the same level of approximation 
as Ql. Both L2 and NLG were fit to describe different 
nuclear properties than those corresponding to the sets 
Wl, Ql, Q2, Gl and G2 ||. They also differ from these 
in the number of nuclei used to fit the parameters. The 
values of the constants used in this paper are given in 
Table ^, and the fitting procedure is discussed in more 
detail later in this paper. 

Other approaches have been applied to studies of nu- 
clei far from stability with various degrees of success. 
Neutron and proton skin formation, the reduction of the 
spin-orbit interaction and the realization of pseudo-spin 
symmetry are some of the phenomena studied recently 
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by various authors 17, Shell-model calculations 
using the CD-Bonn effective potential modified by the G- 
matrix and renormalized using the Q-box method have 
been used to study the importance of the ^So and 
partial waves in reproducing the constancy of the 2^—0^ 
level spacing which is apparent in Sb nuclei ||l9| . Some of 
the methods that have been tested on these nuclei with 
various degrees of success are the density dependent rel- 
ativistic Hartree (DDRH) [p^ |20| , relativistic continuum 
Hartree Bogoliubov (RCHB) ||17|, relativistic Hartree Bo- 
goliubov (RHB) Hartree-Fock-Bogoliubov PH] and 
relativistic Dirac-Hartree-Bogoliubov (DHB) |22| |. In 
these methods pairing correlations have been taken into 
account using a phenomenological force like the Gogny 
force jlHl [itI or a density-dependent zero-range inter- 
action PQI . Other authors [£0|, ^ used the BCS ap- 
proach with a constant pairing strength to take into ac- 
count the effects of the continuum. Perhaps the most 
self-consistent approach has been the use of the RHB, 
in which ph-interactions and pairing correlations have 
been taken into account in a unified way by solving self- 
consistently the RHB-equations with the Klein-Gordon 
equations that describe the meson fields jl^, Other 
authors have derived new parameter sets of known po- 
tentials, e.g. Nilsson [^, Groningen ||2^, to fit the new 
experimental data of binding energies, single-particle lev- 
els, and nuclear radii for these nuclei far from stability. 
Relativistic mean field models (RMF), which have been 
successful in describing bulk properties of stable nuclei, 
are being tested in these new regimes. Some of the mod- 
els use the a — u lagrangian with linear and non-linear 
a interactions Applications of DFT to the nu- 

clear many-body problem have been carried out using 
the QHD-I and QHD-II lagrangians H, |||. Here the 
exchange-correlations are taken directly from the applica- 
tions of many-body techniques in infinite nuclear matter, 
but extensions to nuclei far from stability have not yet 
been performed in this approach. Direct application of 
the QHD-H lagrangian to the doubly-magic Sn isotopes 
to calclulate single-particle level spectrum and nucleon 
densities can be found in . 

The combination of the fundamental idea of DFT and 
the realization that an approximation to the exact en- 
ergy functional can be derived using EFT methods is a 
new approach to the nuclear many-body problem. As 
such, it is still in development. To the extent that a suf- 
ficient number of terms are retained in the effective la- 
grangian to provide the correct energy functional, one has 
actually developed the effective field theory for QCD at 
the nuclear density domain. Extensions to other density 
regimes still has to be explored. This paper explores how 
additional terms in the effective lagrangian contribute to 
the results for ground-state properties of selected nuclei 
far from stability. The EFT methods in this chiral ef- 
fective lagrangian approach have already shown that the 
coupling constants are all natural [Q, therefore the rel- 
evant degrees of freedom have been correctly taken into 
account and the NDA works. Knowing how well this 



EFT approach behaves at these extreme cases of nuclei 
far from stability will give a firmer ground for the the- 
ory. Predictions for the properties of these nuclei is also, 
of course, of direct relevance to experiments, both ongo- 
ing and planned [ p8| , and to element formation in super- 
novas. 

In this paper we see that the increasing level of ap- 
proximation approaches the total binding energy of the 
selected nuclei. There is also very good agreement with 
the binding energies of neighboring nuclei, which differ 
by one particle or one hole from the doubly-magic nu- 
cleus, i.e. the chemical potential. The spins and parities 
of these nuclei, to the extent to which they have been 
measured, are also correctly given. In addition, we show 
the calculated single-particle excitation spectra and ex- 
amine the effects of additional levels of approximation. 
Finally, we present the nucleon densities for neutron rich 
and neutron deficient nuclei predicted with this approach. 

The following two sections explain the effective la- 
grangian used in this paper and describe the calculations, 
and in the last two sections results are presented and con- 
clusions are drawn. 



II. CALCULATIONS 

The chiral effective lagrangian used here was con- 
structed using the relevant lowest lying hadronic degrees 
of freedom The pion field is introduced as the phase 
of a chiral rotation of the identity matrix in isospin space. 
A field ^{x) is defined in the following way: 

^{x) = exp{i7r{x)/f^). 

where the pion field Tr{x) is defined by tt{x) = • tt 

Nucleons are included as isospinor fields, N{x), in 
which the upper component corresponds to the proton 
and the lower component the neutron. 



N 



p{x) 
n{x) 



To include a contribution to the symmetry energy of nu- 
clear matter, a p^-field is added to the lagrangian. 

A non-linear realization of chiral symmetry is consid- 
ered following the work of Callan, Coleman, Wess, and 
Zumino A description of the incorporation of chi- 

ral symmetry in the nuclear many-body lagrangian can 
be found in j^, 0|. In this formalism, the non-linear re- 
alization of the chiral symmetry is defined by a global 
transformation, L and R, of the subgroups SU{2)l and 
SU{2)r such that 



where 



^'(.t) = L^{x)h\x) = h{x)^{x)R^ , 
P'f,{x) = h{x)p^{x)h^{x) , 
N'{x) = h{x)N{x) . 



(1) 

(2) 
(3) 
(4) 
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here p^{x) = ■ Pi_,. 

In addition to the above fields, an electromagnetic field 
Af^ is included and the electromagnetic structure of the 
nucleons has been taken into account. Their anomalous 
magnetic moments, labeled as Ap and A„ for proton and 
neutron respectively, then enter. An isoscalar- vector field 
VJj is needed to correctly represent the bulk properties of 
nuclei. This field is treated as a chiral singlet and its mass 
is taken to be equal to the w-meson mass. Finally, an 
isoscalar-scalar field cj), again a chiral singlet, is included 
to incorporate the observed mid-range N-N attraction. 
Its mass is determined by fitting the calculations to ex- 
perimental results. 

The various interaction terms present in the effective 
lagrangian have been arranged following the NDA from 
Georgi and Manohar |l^ . Details concerning the con- 
struction of the chiral effective lagrangian can be found 
in j|] where the full effective lagrangian is developed. 

In order to solve the equations of motion derived from 
this lagrangian we apply the relativistic mean-field ap- 
proximation in which the meson fields become expecta- 
tion values and thus are treated as classical fields. In 
this approximation the pseudo-scalar pion field does not 
contribute because it has no expectation value. The in- 
terplay between these fields, particularly the scalar and 
vector fields, accounts for the bulk properties of nuclei as 
is well described in Q . For spherically symmetric sys- 
tems, as is the case for double-magic nuclei, the spatial 
components of the vector field vanish. The rho isovec- 
tor field develops only the non-charged component, po, 
since conservation of charge is imposed. The resulting 
equations of motion consist of a Dirac equation for the 
nucleons and a system of non-linear, coupled differential 
equations for the meson fields. 

The Dirac hamiltonian for the nucleon fields takes the 
form: 



h{x) 



-id-V + W{x) + -T3R{x) + f3{M - $(2;)) 
+ l{l + rs)A{x) 



I 

2M 
1 

2M2 



Pa-ifp-TsVR + f^VW) 



{l3s+PvT3)V^A 



2M ^ 



(5) 



Here A = 5Ap(l + T3) + iA„(l — T3) and the numeri- 
cal values used for the anomalous magnetic moments are 



A„ 



1.793, A„ = -1.913. 



The mean meson fields are denoted by W — gvVo, 
$ — 5s 00, R — gpbo, and A = cAq respectively. The 
quantities fp,fv,l3s,l3v are parameters fit to experiment 
(see below). The meson field equations become the fol- 
lowing: 



K3 



M ^ 2 



K4 <I> 



^ .ml. 



,9s/ I \ '"-u Tjr2 



•y p 



2Mg\ 



(6) 



Here ps is the baryon Lorenz scalar density and gs, m^, 
K3, K4, Tyi, rj2, gv, gp, cn, 0.2 are again parameters fit to 
experiment. Furthermore, 

- V^W -h mlW = gl[pB(x) + ^V.(pi(x)f)] 

^ 2 m' M " 

3! ' 



e^gv 
3g^ 



Pchgix) 



(7) 



Here is the baryon density, p]^ is the baryon tensor 
density, pctg is the charge density and Co are param- 
eters. In addition, 
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VpijmiR- ^-^Pchgix) (8) 



M 



9i 



Here pa and pj are the isovector densities, /p, r/p, gp are 
parameters and (7-y=5.01, is the coupling of the photon 
to the w-meson. Finally, 



V^A = eVchg(a;) 



(9) 



where pchg is the charge density. To solve the Dirac equa- 
tion and obtain the single-particle energy levels one has 



hlpa{x) = Ea1pa{x) 

The Dirac solutions are expressed as. 



' (^) ^ — Kr) 



(10) 



(11) 



where the $Km is a spin spherical harmonic and t is equal 
to 1/2 for protons and -1/2 for neutrons. Inserting this 
into Eq.([l0|) and using Eq.(||) a set of two coupled first- 
order differential equations for the G and F functions is 
obtained: 



^ + -) Gair)-[Ea-U,{r)+U2{r)]Fa{r)-UsGa{r) = 
ar r / 

(12) 

^--] Fa{r)+[Ea-Ui(r)-U2{r)]Ga{r)+U3Fa{r) - 
ar r J 



(13) 
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where the single-particle potentials are given by 



Ui{r) 
U2{r) 



W{r)+taR{r) + {ta + ]^)A{r) 



1 



:{(is + 2ta[3v)V^A{r) 



(14) 
(15) 



M - $(r) 

+ A'{r)[{\p + A„)/2 + ta[\p - A„)]} (16) 



:{fvW'(r)+tafpR'{r) 



Here the prime indicates a radial derivative; e.g. W (r) = 
dW{r)/dr. 

Once the single-particle wave functions are calculated, 
the various densities that appear on the r.h.s. of the 
meson equations can be obtained. They are defined as 
follows: 



Ps{x 
Pb{x 
Pb{x 



E 



47^7-2 



{Gl{r)-F^{r)) 



E^2G.(r)^.(r) 

a 

occ ^ . 

E^(2ia)(G^M+^^a^M) (20) 



(17) 
(18) 
(19) 



E 



2ja + 1 



{2ta)2G air) Fair) 



(21) 



where the sum goes over the occupied orbitals. In ad- 
dition to these, the charge density is composed of two 
pieces: a direct nucleon charge density and the vector 
meson contribution 

Pchg = P<i(x) -I- Pm(x) (22) 

The direct part is given by 

Pdix) = ppix) + ^V-iplix)r) 

+ ^i/^s^^PBi^) + /3.V2p3(x)] (23) 

and the vector meson contribution arising from the cou- 
pling of the neutral vector mesons to the photon (VMD) 
takes the form 



Pmix) 



1 



1 



(24) 



Here the point proton density pp and nucleon tensor 
density pj are given by 



1/ 

Pp 2^PB + P3) 



~ ipl.ix)iXpd ■ ripa-i^ 



(25) 
(26) 



To solve the meson equations a Green's function 
method is used. The equations are iterated until a consis- 
tent solution for all the meson fields is obtained, that is, 
we obtain a O^'^-order solution to the meson equations by 
neglecting all non-linear terms on the r.h.s., and we then 
use this to obtain a I'^^-order solution, continuing until 
convergence is achieved. Because NDA guarantees that 
each additional term on the rhs of Eqs. (^ - (^ is smaller 
than the previous one, convergence is both expected and 
obtained. 

The whole system of Eqs. (|l|), (|l|) and (|) - (|) 
is solved self-consistently until a global convergence is 
reached. 

The values of the various constants used in this paper 
are given in Table § with a label indicating the parameter 
set. These values are taken from |Q, ^. The constants 
were obtained by fitting several nuclear properties [0. 
From the 16 constants in Q the last 3 where determined 
by fitting to electromagnetic properites of the nucleon. 
The sets labeled Wl, Ql, Q2, Gl and G2 were obtained 
using a weighted generalized defined by 



-EEi 

i X 



.X, 



cxp 



X. 



(i) 

th l2 



to fit a total of 29 observables listed as follows: 

• The binding energies per nucleon E/B with 
W=0.15%; 

• The rms charge radii {r'^)^i^^ with W—0.2%; 

• the diffraction-minimum-sharp, d.m.s., radii Rdms 
with W=0.15%; 

• The spin-orbit splittings AEso of the least-bound 
proton and neutrons with W=5% for -^^O, 15% for 
208pb, 25% for 40Ca and ^^Ca, and 50% for ^^Sr; 

n 

• The proton energy EpilhQ/2) ^-nd the proton level 
splitting Epi2d3/2) - Epilhii/2) in ^ospb with 
W=5% and 25%, respectively; 

• The surface-energy and symmetry-energy deviation 
coefficients 6a2 and 5a4 with a weight W=0.08. 

The parameter set L2 was obtained by fitting the con- 
stants to reproduce the infinite nuclear matter proper- 
tics @: 

• Saturation density p°g = 0.1484 fm"^. 

• Binding energy per nucleon E/B = -15.75 MeV. 

• Bulk symmetry energy = 35 MeV. 

• rms charge radius of ^°Ca — 3.482 fm. 
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TABLE I: Values of the different parameters in each parameter set used to solve the self-consistent mean-field equations 
from |. rus is in MeV. 



Constant 


L2 


NLC 


Wl 


Ql 


Q2 


Gl 


G2 


rUs 


520 


500.8 


566.26 


504.6 


509.6 


506.7 


520.3 


gl 


109.63 


95.11 


138.93 


103.6686 


97.7699 


97.39 


110.16 




190.43 


148.93 


203.97 


164.6963 


149.2 


147.09 


162.88 


up 


65.23 


74.99 


95.551 


77.9558 


73.2256 


77.033 


89.936 


m 












0.07060 


0.64992 














-0.96161 


0.10975 


K3 




1.9194 




1.6582 


1.7424 


2.2067 


3.2467 


Ki 




-7.3923 




-6.6045 


-8.4836 


-10.090 


0.63152 


Co 










-1.7750 


3.5249 


2.6416 














-0.2722 


0.3901 














1.8549 


1.7234 


a2 












1.788 


-1.5798 


fv 












0.4316 


0.6936 


fp 






3.7328 


4.1328 


4.2640 


4.1572 


3.8476 


Ps 






-0.38482 


-0.10689 


0.01181 


0.02844 


-0.09328 


13^ 






-0.54618 


-0.26545 


-0.18470 


-0.24992 


-0.45964 
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FIG. 1: Total binding energy of the selected nuclei. The left 
hand side shows two models with the same level of approx- 
imation but fit to different data sets (see text). The right 
hand side shows models with increasing level of approxima- 
tion. NLC and Ql have the same level of approximation as 
do Gl and G2. 



III. RESULTS 



In this section we show the results of the calculations 
of several ground-state properties of the nuclei ^|QSn82, 
^5oSn5o, ||Ni2o, and 28^150 and of the neighboring nuclei 
differing by one particle or one hole. We start by con- 
sidering results for the total binding energy, which are 
shown in Fig. |l| 

The plot is divided in two parts which show differ- 
ent forms of the convergence of the models. On the left 
are the results of total binding energy using the simplest 
lagrangian, i.e. with no non-linear terms. The two pa- 
rameter sets L2 and Wl represent the result of fitting 



the constants of the model to different nuclear proper- 
ties. As described in the previous section, the Wl set 
is fit to more data than L2. The improvement of the 
results is clear. The experimental result is shown as a 
long-dash line in the plot. The error bars are of the same 
size as the thickness of the line showing the experimen- 
tal value. The experimental results are taken from [30|, 
except for 28^^120 where we use the value quoted in [31|. 
On the right side of Fig. are the results obtained from 
the more elaborate lagrangians, which in general include 
different non-linear interactions as shown in the previous 
section (see Eqns. (||) - (^). The first two points, corre- 
sponding to the NLC and Ql parameter sets, represent 
the same level of approximation of the lagrangian. We 
see here, again, that the use of an extended set of observ- 
ables in the fitting procedure dramatically improves the 
results. In this case, the Ql set already gives an excellent 
result for all nuclei. 

As we move to the right on this part of the plot, we 
see the effect of additional interaction terms in the la- 
grangian. The inclusion of a fourth-order vector field 
term in the lagrangian, see Eq.(|^) and Table |, gives a 
larger binding energy. This effect is further balanced by 
additional non-linear terms in the Gl and G2 sets. At 
this level, the sets Gl and G2 give almost the same results 
and cannot be differentiated. Since both sets correspond 
to the same level of approximation, we are looking here 
at the sensitivity of the results to the variations of the 
fitted constants. The same features are verified for all 
the nuclei tested. Figure ^ shows the percentage devia- 
tion in estimating the total binding energy of the selected 
nuclei. We observe that it lies below approximately 1% 
in all cases. In all cases the parameter set Ql gives the 
smallest deviation. Extending the parameter set to Gl 
and G2, while validating the NDA and the assumption 
of " naturalness" as discussed in [1] , does not significantly 
improve the quality of the fit to this quantity for the nu- 
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clei shown. For the four doubly-magic nuclei tested we 
obtain better agreement for the case of the Sn-isotopes 
than for those of Ni. 

For the Sn isotopes, we extended the calculations of 
the total binding energy to cover the entire range of even- 
even nuclei. Figure ^ shows the results for the G2 param- 
eter set for the isotopes ^gpSnso - -^|QSn82. The percent- 
age deviation in this extended case is shown in Figure ^. 
Here we observe that the deviation for the three sets Ql, 
Gl and G2 is in absolute value less than 1%. In general, 
all three parameter sets give less binding energy than 
observed. For all parameter sets shown in this figure, 
Gl gives overall the smallest deviation, except for the 
two nuclei ^jpSngo and ^5QSn82, which lie at the extreme 
edges of the plot. This would imply that there is indeed 
improved convergence overall in going from Ql to the ex- 
tended parameter sets Gl and G2, and Gl gives better 
results than G2 over the entire range. 

From this figure we can observe the effects of the ac- 
curacy of the isovector contributions to the effective la- 
grangian. Following, for example, the Gl results we ob- 
serve that in the center of the plot there is essentially a 
constant deviation from the experimental values, while 
this increases close to the edges of the isotope range. A 
study of the problem of the isovector contribution to the 
description of finite nuclei has been presented recently by 
Furnstahl [Ei. 
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FIG. 2: Percentage deviation of total binding energy of the 
four selected doubly-magic nuclei. Parameter sets Ql, Gl and 
G2 are shown. 

We can conclude that a more sophisticated lagrangian, 
fit to stable nuclear properties, can reproduce ground- 
state total binding energies of nuclei far from stability. 

It is important to note that even though the constants 
of the model were fit using only stable doubly-magic nu- 
clei, the results for binding energies of these nuclei far 
from stability are remarkably good. 

Figures || through |l0| show the level structure of neigh- 
boring nuclei differing by a single particle or hole, relative 
to the doubly-magic nuclei. They are separated in iso- 
topes and isotones. The core nucleus is plotted in the 
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FIG. 3: Comparison between experimental and calculated to- 
tal binding energies for Sn-isotopes using the G2 parameter 
set. 
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FIG. 4: Percentage deviation of the total binding energy for 
Sn-isotopes using Ql, Gl and G2 parameter sets. The stable 
isotopes are indicated in the plot. 



center of the plot marked with 0+ . On the left we have 
put neighboring nuclei with one more particle. On the 
right the corresponding nucleus with one less particle is 
shown. When available, we include the experimental val- 
ues. These results were obtained using the G2 param- 
eter set, although the Gl set gives very similar ground 
state binding energies. In all these cases there is excellent 
agreement with the experimental data pO| . 

In the case of ^lo^ngi the above reference gives as 

ground-state a value of | , but recent experimental re- 
sults from Bhattacharyya, et al. |^ indicate that the 
ground-state is in fact ^■ in accord with our results. 

We have applied the chiral effective lagrangian to cal- 
culate the ground-state properties of selected doubly- 
magic nuclei. The inclusion of additional interaction 
terms in the effective lagrangian, consistent with NDA, 
gives results that approach asymptotically to the exper- 
imental value for the total binding energy. We extended 
this calculations to neighboring nuclei, differing by one 
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particle or one hole, and determined the binding energy 
relative to the doubly-magic one and found an overall 
agreement of better than 1% for the total binding energy 
of the doubly-magic nuclei and better than 10% for the 
relative binding energy of the neighboring nuclei. In this 
latter case, the spin and parities of the ground-states are 
correctly predicted. 

Single-particle excitation spectra for the nuclei ^InSnga 
and 1^52 Sbg2 have been measured by Hoff, et al. and 
by Sanchez- Vega et al. |3j]. The level spectra for these 
nuclei is shown as a function of parameter set in Figs. |ll| 
and |l2|. Here we chose only the G2, Gl and Ql sets since 
these give the best results for the total binding energy of 
i^|QSn82. From the discussion of the DFT approach we 
do not expect these single-particle and single-hole exci- 
tation spectra to be well reproduced, though we observe 
that there is a slow convergence of the results with higher 
levels of approximation. This might indicate the need for 
many more terms in the lagrangian. On the other hand 
the excellent agreement shown in Fig ^ for the total bind- 
ing energy, and Figs. |^ and ^ for the ground-state prop- 
erties (i.e. chemical potential leading to the neighboring 
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FIG. 12: Comparison between single-particle level spectrum 
of ^IgSnga with recent experimental data 
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single-particle and singlc-liolc nuclei) suggests that the 
current level of approximation is enough to reproduce 
the ground-state observables the theory is designed for. 
We observe that increasing levels of approximation give 
better results in the level spacing, e.g. Fig |ll| though 
the ordering is essentially independent of the level of ap- 
proximation. Note that in the case of Sb82 we do not 
include results for the Ql parameter set for levels above 
2d5/2. This is because the program used did not find 
those states to be bound, a problem that is not present 
in the case of ^IgSuga. This feature is expected since 
the higher order terms included in the Gl and G2 pa- 
rameter sets do not only contribute to the bulk proper- 
ties of the nuclei but also improve the description of the 
single-particle spectra. In Figs. |l^ and |l^ we show the 
results for the neutron and proton densities of ^|oSn82 
and ^5oSn5o. In both cases the density of neutrons is 
always higher in the interior of the nuclei. This is par- 
ticularly interesting in the case of ^^qSusq since there are 
an equal number of neutrons and protons implying that 
the protons are pushed out, due to the Coulomb poten- 



tial, forming a thin layer at the surface. In contrast, in 
the neutron rich case the neutrons extend far out and 
form a neutron layer. In Fig. |l5| we compare the neutron 
densities of these two doubly-magic Sn nuclei. These re- 
sults are in good agreement with calculations made using 
other approaches p5[ |. 

The calculations of the nuclear properties using the 
DFT approach do not include the contributions of pair- 
ing interactions. Even though most modern calculations 
use some type of pairing, as mentioned before, the goal of 
this study is to determine the convergence of the results 
for the optimal observables as new interaction terms are 
included. Extending these calculations beyond these ob- 
servables gives less satisfactory results, as has been seen 
in the single-particle and single-hole excitation spectra of 
neighboring nuclei, Figs[l| and |l^. Although beyond the 
original purpose of the paper, it is interesting to compare 
the experimental values of the two- neutron separation en- 
ergy S2N for the series of nuclei ^^Snso - ^fpSngg with 
our calculated results. We give results for G2 only, since 
Gl gives essentially the same values. To make a rea- 



10 




FIG. 14: Proton and neutron densities for sgSnso 
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TABLE 11: Results of extracting the configuration energy 
from the experimental data on S2N using the pairing model 
described in the paper. 
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methods. In one case we put only one particle in the 
j-th sub-shell, while the other lower shells and sub-shells 
where filled, and determined the total binding. Extract- 
ing the contribution from the lower filled shells and sub- 
shells, to the total binding energy, gave the configuration 
energy. In the second case the j-th sub-shell was com- 
pletely filled. We proceeded as before and caclulated the 
total binding energy and subtracted from it the contribu- 
tion from all the lower filled shells and sub-shells. This 
final result was averaged by dividing it by the total num- 
ber of neutrons in that j-th sub-shell, 2j + 1. Both meth- 
ods gave essentially the same result. Figure |l^ shows the 
comparison between the extracted and calculated con- 
figuration energies. The plot shows that, though not 
quantitatively, the present approach reproduces the level 
ordering. 

The inclusion of pairing correlations is a subject that 
goes beyond the scope of the present study, but it is rec- 
ognized that some sort of pairing has to be included to 
obtain a more complete description of nuclear observ- 
ables. Work on this subject has shown that in the case 
of standard RMF models the approach using the DHB 
seems to be the most straight forward way to include 
such effects The role of pairing in this EFT ap- 

proach is currently under investigation by the author. 
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IV. CONCLUSIONS 



FIG. 15: Comparison between neutron densities of goSng2 
and ^soSnso. 

sonable comparison, we extract the configuration energy 
contribution coming from the mean-fields by eliminating 
from the experimental data the pairing contributions. To 
do this we use a simple solvable pairing model [ p5| of j- 
dependent strength Gj corresponding to a pure pairing 
force. The result of the model is that S2N is related to the 
number of particles in the j-th sub-shell by the formula: 

S2N = 2egff(j) + Gj{2j + 5-2n) 

where n is the number of particles in the j-th shell. Fit- 
ting the experimental data to this formula, which indeed 
goes through all the experimental points, gives the val- 
ues of the configuration energy 2eg-g-(j) and Gj shown in 
Table 0. 

To obtain the calculated configuration energy Seff{j) 
using the present approach, we calculated it using two 



From the results shown in this paper we draw the fol- 
lowing conclusions: 

1. As additional interaction terms are included in 
the effective lagrangian, following the prescriptions 
given by naive dimensional analysis, the calculated 
results of ground-state properties of selected nuclei 
far from stability converge to the observed exper- 
imental values. The overall agreement, averaged 
over all tested nuclei is better than 1% in the total 
binding energy. The relative ground-state binding 
energy of neighboring nuclei, differing by one parti- 
cle or one hole, are well reproduced, with a relative 
error of better than 10%. The spin and parity of 
these nuclei, where measured, is also correctly re- 
produced. 

2. The effect of fitting the constants of the Chiral Ef- 
fective Lagrangian to a large number of observables 
greatly improves the results for the total binding 
energy of these nuclei far from stability. 
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3. Nucleon densities calculated with this approach are 
in agreement with previous works. 



4. The effective field theory and density functional 
theory approach provides a predictive framework 
for dealing with QCD in the strong-coupling, nu- 
clear physics regime fl^ . 



5. Additional interactions to reproduce the effect of 
pairing are needed to extend the theory to repro- 
duce additional ground-state properties. 



FIG. 16: Comparison between configuration energies for Sn- 
isotopes. The dashed line corresponds to calculated values 
using the G2 parameter set and the solid one is the value ex- 
tracted from the experimental data on two-neutron separation 
energies 5'2jv using the simple pairing force model described 
in the paper. 
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